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HYDRODYNAMIC IMPACT OF A SYSTEM WITH A SINGLE ELASTIC MODE 
I—THEORY AND GENERALIZED SOLUTION WITH AN APPLICATION 
TO AN ELASTIC AIRFRAME? 


By WILBUR L. Maro 


SUMMARY 


Solutions of impact of a rigid prismatic float connected by a 


massless spring io a rigid upper mass are presented. The 


solutions are based on hydrodynamic theory which has been 
ezperimenlally confirmed for a rigid structure. 

Equations are giren for defining the spring constant and the 
ratio of the sprung mass io the lower mass so that the two-mass 
xystem provides representation of the fundamental mode of an 
airplane wing. The forces calculated are more accurate than 
the forces which would be predicted for a rigid airframe since 
the effect of the fundamental mode on the hydrodynamic force 
ig taken into account. The response of the two-mass system 
gites the response of the represented mode and, although no 
proviston te made for taking into account the effect of secondary 
modes on the hydrodynamic force, means are indicated whereby 
the results may be used to approximate the response of modes 
other than the fundamental mode. 

Time histories of the hydrodynamic force and structural 
response are giren for wide ranges of mass disiribulion and 
ratio of natural period to the period of the impact. By use of 
nondimenstonal coefficients these results are made applicable 
fy different combinations of relocity, weight, angle of dead rise, 
and fluid density. Although the equations permit solutions for 
different combinations of flight-path angle and trim, an approxi- 
mation ts giren for correcting the results for the combination for 
which solutions are diren to other conditions unthin a narrow 
range indicated to be of primary interest to the design engineer. 

In a comparison of the theoretical data with data for a serere 
flight-test landing impact, the effect of the fundamental mode on 
the hydrodynamic force is considered and response data are 
compared with experimental data. Consideration of the funda- 
mental mode alone fails to account for the fact that during the 
impact partial failure of the inbuard-engine mounts occurred, 
but use of the theoretical solutions to approximate the effects of 
further wing torsion leads to substantial agreement. 


INTRODUCTION 


In recent years the development of lerge airplanes has 
caused the elastic behavior of airframe structures during 
landing impact to become important. The work which has 
been done on this problem has been handicapped by lack of 
proper knowledge of the time history of applied ground 
reaction. This situation has been particularly acute for 


seaplanes because of difficulties. in measuring the hydro- 
dynamic force, the seaway, and the manner of contact with 
the sea way. 


In order to facilitate the interpretation of flight data and — 


to lead to the prediction of design loads on & rational basis, 
a theoretical hydrodynamic study was made and tests of a 
rigid float were conducted at the Langley impact basin. 
Since the results of these tests agree with the theoretical 
results for wide ranges of the pertinent variables in numerous 
force time histories (reference 1), it is &ssumed that the 


theory may also be used in considering the effect of the 


upper-structure elasticity of a seaplane on the motion and 
force characteristics of the hull proper, which is assumed to 
be rigid. 

The bending of wings during impact, which for modern 
flying boats is the primary structural action, is considered in 


the present report by reducing the fundamental mode to an 


equivalent two-mass system. The results are presented in 
& form suited to general application and are compared with 
experimental results for a particular case. The equations 
showing the method of solution are included in appendix A, 
and a sample data sheet is given as table I. 


SYMBOLS 

ti time required for one-fourth cycle of natural 
vibration 

T time between initial contact and maximum hydro- 
dynamic force for rigid structure 

t time elapsed after initial contact 

Mr lower, or hull, mass of two-mass system 

Ms upper, or sprung, mass of two-mass system 

m gross mass (1V/[g or ms--m;) 


W gross weight 

g acceleration due.to gravity 

K spring constant of spring connecting ms and mz, 
force per unit deflection 

n, '  &cceleration normal to water surface of nodal point 
of elastic system, multiples of acceleration of 
gravity; for two-mass system, acceleration of 
center of gravity 

fio oscillatory acceleration of hull about center of gravity 
of two-mass system or nodal point of represented 
mode, multiples of the acceleration of gravity - 


t Supersedes NACA TN 1398, “Solutions for Hydrodynamic Impact Force and Response of a Two-Mass System Nun an Application to an Elastic Airframe” by WHbur L. Mago, 1947. 
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Vo resultant velocity at instant of contact. with water 
surface 

p mass density of fluid 

T angle of trim; angle of hull keel with respect to plane 


of water surface 

Yo flight-path angle at contact; angle between flight 
path and plane of water surface 

B angle of dead rise 


E 
C, nondimensional time coefficient (: Vo (88) ) 


g? 2A 


1/8 
C, nondimensional draft coefficient CM (sd ) 


Yamar draft at instant of maximum acceleration . 
DP n&tural bending frequency 





C; nondimensional load-factor coefficient (5 =E, 


Where units are not given, any consistent system of units 
may be used. 


THEORY 


HYDRODYNAMIC 


The hydrodynamic theory used in the present report is the 
same as that developed in references 1 and 2.. A basic 


differential equation which gives the instantaneous force in. 


terms of the instantaneous position and motion of the float 
is given in reference 2. This equation is used herein to 
determine the effect of airframe elasticity in altering the 
motion and force time history (appendix A). The solution 
is based on the assumption that the float does not change 
trim during impact. In this connection the pitching moment 
may be large, but the time of the impact is short enough to 
warrant (at the present stage) neglect of the resulting angular 
velocities and displacements. | 

The solution presented herein is | for a prismatic float 
with such beam loading that the chines do not immerse, 
during impact. For waves that give the severe design 
condition of full-length impact, conventional beam loadings 
are small enough to cause the maximum force to occur at 
drafts sufficiently small to make the effects of finite width 
and chine flare secondary. Reference 1 indicates that for a 
conventional float neglect of the pulled-up bow is justified 
when the trim is 3° or greater. Although for high-trim 
landings initial contact by the afterbody may substantially 
change the trim before (hae. main forebody impact, the 
neglect of afterbody loads is justified because, during the 
main impact, the shielding of the afterbody by the forebody 
due to depth of the step and to keel angle is such as to mini- 
mize the importance of afterbody loads. 


STRUCTURAL 


A simplified representation of primary elasticity of an 
airframe is shown in figure 1. A rigid lower mass my is 
considered to be connected by a massless spring to a rigid 
upper mass ms. In determining the fundamental bending of 
airplane wings part of the wing mass must be included in 
mz, and part of the wing lift should be applied to mz. In 


the present report the gravity force on each mass is assumed | 


to be balanced by wing lift. 
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FiGUuRE 1.—81!mplifled representation of primary elasticity of an alrframó. 


The problem of determining the properties of the two- 


mass system so that it is representative of the primary 


elastic action of the airplane is rather simple if it is assumed 
that during the impact the structure deflects with the shape 
of its fundamental mode of vibration. The requirements are: 

(1) The total mass of the simplified system must equal 
the total mass of the airplane in order that the proper nodal 
or center-of-gravity accelerations can be obtained. 

(2) The energy of vibration for the same amplitude of 
the hull and lower mass (relative to the nodal point) must 


be the same for the two-mass system as for the considered. 


mode of the airplane structure. 

(3) The natural frequency of the two-mass system must, 
be the same as the frequency of the considered mode of 
airplane vibration. 

Equations which permit determination of the masses and 
spring constant of the simplified system so that it meets 
these requirements are given in appendix B. These equa- 
tions and the foregoing requirements are applicable for both 


 lendplanes and seaplanes. ' 


In the present report the represented structural mode is 
considered to be devoid of vibration prior to the instant of 
impact. Thus, the computations may represent either & 
first impact or a subsequent impact resulting from a bounce 


. sufficiently high to cause aerodynamic and structural 


damping to stop the vibration during the time the seaplane 
is in the air. This report does not give & representation of 
successive impacts, such as might occur in seaway, which 
lead to accumulative or resonant effects. Available flight 
data indicate that a single heavy impact, such as that 
considered herein, is the primary cause of structural failures. 


The response of the two-mass. system is obtained in eon- | . 


nection with the calculation of the time history of the hydro- 
dynamic force, and from this result the complete response 
of the. represented mode can be obtained by the simple 
procedure ‘given in appendix B and demonstrated in the 
section entitled "Comparison with Experiment." The 
response of other modes to the force computed on the basis 
of the fundamental mode can be separately determined 
and superposed (reference 3). In order to minimize the 
complexity of the solution, however, the present investi- 
gation does not provide for taking into account the effect 
of the other modes on the hydrodynamic force. Although 


the other modes may have a substantial effect ọn the local | .. 


loads in the structure, the effect of these modes on the 
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hydrodynamic force is considered to be secondary as com- 
pared with the effect of the fundamental mode. 

If a large number of solutions for the two-mass system 
have been made in order to determine the effect of the funda- 
mental mode of different wings on the hydrodynamic force, 
the response of modes other than the fundamental can be 
upproximated from the use of a solution for a mass ratio 
and ratio of the time period of the impact force to the natural 
period of the two-mass system representative of the con- 
sidered mode. An example of such use to approximate the 
effects of wing torsion is given in the present investigation 
. in a comparison of results of computations with experi- 
mental results. 

If the response of more than one mode is considered, the 
structural and aerodynamic damping, which are not con- 
sidered herein, are important factors in determining the ex- 
tent to which the maximum response of the different modes 
should be superposed without regard to phase relationship. 
It is expected that the effect of the damping will be most im- 
portant for the higher modes and that a result leading to con- 
servative design will be obtained if damping is not considered 
and the maximums of the first two or three modes are super- 
posed without regard to phase relationship. | 


RESULTS 


Solutions of the equations in appendix A were made for 
wide ranges of the pertinent variables. Time histories of 
the calculated nodal acceleration. or hydrodynamie force in 
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Fracre 2.— Variation of force time history with time ratio. 
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terms of the weight, are given in figures 2, 3, 4, and 5 for. 


ratios of the sprung mass to the lower mass equal to 0.25, 
0.60, 1.00, and 1.36, respectively. Each figure is three- 
dimensional; the third dimension is /,/f; which is & ratio of 
the period of natural vibration to the speed of the impaet. 
In representing the period of natural vibration, £, is taken as 
the time required for one-fourth of & cycle. "The speed of 
the impact is represented by making t; equal to the time be- 
tween initial contact and maximum acceleration for a rigid 
structure. If the time to reach maximum force for the elastic 
structure should be used in defining ts, discontinuities in the 
time to reach maximum force would cause discontinuities in 
the time-ratio scales of the plots. (See figs. 2 to 5.) 

An expréssion for t£, may be obtained from the relation 


LL f, and equation (A4) in appendix A. The equation 


for t: is as follows: 


t Q (2 HS l ` 1 
L a pg ( ) 
where l l 
Ca Ume coefficient at instant of maximum acceleration __ 
for rigid body (0.678 for 8—22.5?, yp=14°, and 
T=3°) 
The expressions for t, and ¢; may be used to determine that 
ta TVo (e) T RET 
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FIGURE $.— Variation of force time history with time ratio, 0.80 B= 2254°: v9 14?; 72^. 


The oscillatory acceleration is given in figures 6 to 9. 
These figures are the same as the figures giving the nodal 
acceleration (figs. 2 to 5), except that the acceleration plotted 
is the difference between the hull acceleration and the nodal 
acceleration., The time history of the hull acceleration can 
be obtained by summing the two plots for a particular mass 
ratio. On the basis of the assumption that during impact 
the airframe structure deflects in a particular mode, a time 
history of the acceleration of any point in the structure is 
obtained from the results for the equivalent two-mass system. 
by the following procedure: — . | 

(1) From the deflection eurve of the represented mode, 
obtain the ratio of the deflection of the point of interest to the 
deflection of the hull. Both deflections are taken relative 
to the nodal point. 


(2) Multiply this ratio by the oscillatory acceleration 


given either by figure 6, 7, 8, or 9 or by interpolation between . 


these figures for the mass ratio of the equivalent two-mass 
system. 


(3) Add result to the nodal acceleration given by figures 
2 to 5. 


Time histories of the acceleration given in figures 2 to 9 
are on a nondimensional basis. The nondimensional eo- 
efficients, which contain velocity, weight, fluid density, and 
acceleration of gravity, were used in reference 1 10.8. com-, 
parison of theoretical data with impact data for a float 
having an angle of dead rise of 228?. 

Application to other angles of dead rise.—The function of 
the angle of dead rise can also be included in. the nondimen- 
sional coefficients, but in the present investigation this : 
function is isolated and treated as a factor for correcting the 
results presented for angle of dead rise of 224? to other deud- 
rise angles. The pertinent relationships between results for 
different angles of dead rise may be expressed as follows: 


1 1/8 l 
=e rai " 
nyo LF) f(A]? (4) 

where 

f(8) function representing variation of virtual mass for 


two-dimensional flow with angle of dead rise 
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FiGUCRE 4.—Varietion of force time history with time ratio. ——m]1.00; B:22214?; 9a 147; 7937. c 


fid} aspect-ratio factor including effect of angle of dead 
rise on aspect ratio; ratio of virtual mass for 
three-dimensional flow to virtual mass for two- 
dimensional flow 

Np acceleration at any point, either oscillatory, nodal, 


or total 

The shape of the force and acceleration curves for a given 
value of t,/f; is independent of angle of dead rise, but the 
effect of angle of dead rise on t, as given by relation (3), 
does enter into the determination of the value of £,/t, for a 
particular solution. After the value of t,/t, for a particular 
solution has been determined, the acceleration and time 
values for an angle of dead rise of 224° are proportioned by 
means of relations (3) and (4) to the corresponding values 
for the angle of dead rise used in determining f,/t, in order to 
obtain the proper acceleration history. , 

Although adequate impact data have not been available 
for checking the theoretical equations for angles of dead rise 
other than 224?, the theory is equally applicable to planing 
floats. Study of planing data has shown that the functions of 


$15488—54———21 


"EL 


angles of dead rise used in equations herein ere approximately — — 
correct for angles of dead rise ranging from 15? to 30°. The | 
functions are: S 


(35-1) (8) 
TAL aa o) 


Until improved functions of angles of dead rise are obtained, 
functions (5) and (6) should be substituted in equations 

(1) and (2) and in relation (4) to correct for angles of dead 
rise within the range from 15? to 30°. Rough approximation - 
can be obtained by use of functions (5) and (6) for angles of 
dead rise greater than 30? but not for angles of dead rise much 
less than 15°. In reference 2 there is a discussion of the _ 
inadequacy of f(A) for aspect ratios which normally occur 
for small dead-rise angles. For very large angles of dead rise 
and moderate velocity the static forces, which are not 
considered in the present Investigation, become of greater 
importance. 
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Approximate correction to other flight paths and trims.— 
The flight-path and trim conditions of primary interest to the 
designer tend to be independent of the variables in seaplane 
design and constitute a limited range within which approx- 
imate correction of results to different angles of flight path 


and trims can be made without necessity for repeated time- . 


history solutions. For any particular combination of 
horizontal speed, rate of descent, and trim, the most severe 
impact load for most of the structure occurs when the sea- 
way is such that the keel contacts a wave slope approximately 
parallel to it. The effective angle of flight path and trim for 
such an impact are defined relative to the inclined wave slope; 
therefore, the trim which gives maximum force is zero. The 


largest flight-path angle, relative to the keel and to the . 


critical wave slope, is also associated with the most severe 
force. The value of the largest flight-path angle is not so 
definite as the critical trim but tends to be independent of 
variations in size and wing loading. The velocity of the wave 
should be considered in determining the contact speed and 
flight-path angle. | 
The equations and method of solution given in appendix 
A permit solution for different flight-path angles and trims; 
however, approximate correction of the results in figures 2 to 
9 to other positive contact angles can be made by assuming 


that the proportionate effect of the structural elasticity on 
the hydrodynamic force is solely dependent on the ratio 
t./t,. Curves given in reference 1 show values of Ci for 


different flight-path angles and trims, which may be sub- 
stituted in equation (2) to obtain the value of t,/t, for differ- 
ent contact angles. In making the approximate correction, 
the solution presented herein for the obtained value of the 
ratio f/f; should be used to approximate the shape of the 
curve giving the desired time history. The load or accelera- 
tion scale should be corrected to the different contact angles 
by multiplying the present result by the ratio, determined 
from curves given in both references 1 and 2, of the load- 
factor coefficient for the different contact condition to the | 
load-factor coefficient for the conditions of 3? trim and 14? 
flight-path angle considered herein. Correction of the time 
scale involves a similar procedure in which the time coef- 
ficient is used rather than the load-factor coefficient. 

Since the force curves for a rigid body are approximately 
the same shape for different angles of flight path and trims 
(reference 1), the approximate method of correcting to dif- 
ferent angles of flight path and trims would be almost correct 
if the structural elasticity did not affect the hydrodynamic- 
force curve. The. percentage change in the force on the 
float due to elasticity is a function of the percentage change 
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in draft caused by the elastic compression for a given center- 
of-gravity position. An indication of the validity of the 
approximate correction is obtained by studying the extent 
to which the ratio of the spring deflection to the draft is 
constant for impacts of the same values of ¢,/f; at different 
angles of flight path and trims. From the expressions for 
talin Cr, Ca, and C, with the spring deflection assumed to be 
proportional to the hydrodynamic force, this ratio may be 
represented by the expression C;C?/C,, m which all values 
are for the instant of maximum acceleration. The variation 
of this expression with flight path for an impact of a rigid 
float at 3° and 12° trim is given in figure 10. Values of C: 
and C. used in obtaining this figure are given in reference 1; 
values of C4 were obtained in conjunction with the data of 
reference 1 but have not been published. 

In the present report the numerical values of C;C?/C, 
have no significance and they are of interest only because of 
the extent to which they are constant. Figure 10 indicates 
that for large flight-path angles and small trims the ratio is 
approximately constant. The deviation from a constant 
value of this ratio is due to plening forces which exist in an 
oblique impact and become more important for low flight- 
path angles and high trims. The conditions of large flight- 
path angle and small trim previously adjudged to be of 
primary interest to the designer constitute the ranges in 
which the deflection ratio is fairly constant and the approxi- 
mate method of correction should give a fair degree of 
accuracy. The present solutions are considered to be for 
conditions suited to correction of the results to other con- 
ditions of greatest practical interest; they represent a moder- 
ately severe combination of flight-path angle, wave slope, 
and trim chosen to facilitate correlation of the theory with 
an impact which resulted in substantial damage to a well- 
instrumented flying boat during flight tests. 

The equations presented herein are not valid for zero 
trim; an assumption that the float is prismatic gives solu- 
tions of infinite wetted length and infinite force for zero 
trim. Solution for 3° trim and a prismatic float is much 
simpler than a correct solution for 0° trim because necessity 
for consideration of bow shape is eliminated. The solution 
for 3° trim may be taken as an approximation of the critical 
design load or, as illustrated in the following section in & 
comparison of theory with experiment, an empirical factor, 
which includes bow effects, may be used to convert values 
of acceleration and time for 3? trim to values for 0? trim. 


COMPARISON WITH. EXPERIMENT 


Impact basin data have not been obtained for suitable 
models, and most flight landing data have been inadequate 
for the present study. The only data which appeared suit- 
able for this comparison are those which were obtained with 
a four-engine flying boat, the data for which have not been 
published. Data were obtained for a large number of test 
landings, but only one of the impacts is very well suited to 
the present analysis. This impact gave loads sufficiently 
high to cause large effects of elasticity of the wings. The 
impact occurred against the flank of a sizable wave (4 ft) 
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Force-penefratian ratio, 





FIighíi-pa!h angle, Ta, deg 
FiGCRE 10.— Variation of force-penetration pato at maximum force with flight-path angle — 
and trim. 


and thus facilitated the use of results based on a planar water 
surface. The present comparison of theory with experiment 
will be restricted to this impact since other impacts involved 
more complex contact with seaway and gave less force. 

The horizontal speed, rate of descent, and trim were re- 
corded. A large number of pressure instruments distributed 
in the hull permitted determination of the water surface 
relative to the hull. Data recorded by these instruments 
indicated that the wave slope in contact with the hull was 
approximately planar, that the trim relative to the wave 
slope was 0°, that the resultant velocity, considering the speed 
of the wave, was 85 feet per second, and that the flight-path 
angle relative to the wave slope was 14°. 

Structural data available for the test flying boat are not 
adequate for the present analysis. Since this flying boat has 
the same number of engines and approximately the same 
gross weight and horsepower as a landplane for which & 
large amount of structural data is available, assumption is 
made that the flying-boat wing has the same mode shape and 
mass distribution as the wing of this landplane. Use in 
equation (B6) in appendix B of data for the landplane given 
in reference 3 leads to the following mass ratio of the two- 
mass system representing the fundamental mode: 


Ds 0.25 
Mr 





Based on study of the accelerations at several points in the 


' test flying boat during periods of relatively free vibration in 


which the fundamental wing bending mode appeared to be 
predominant, & natural frequency of 3.6 cycles per second 
was selected for use in the present example. The funda- 
mental mode frequency of the wing of the landplane is 3.4 
cycles per second. (See reference 3.) | 
Since the mass ratio is equal to 0.25, figures 2 and 6 are 
used to approximate the action of the fundamental mode. 
Furiher, substitution of values for conditions for this impact 
in equation (2) results in use of the specific time history 
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given for 12, Substitution of the contact conditions in 
d 


the load and time coefficients fixes the Joad and time scales 
in an absolute sense. Correction from 3? trim to 0? trim may 
be made by assuming that the shape of the time history is 
approximately the same for both conditions. Empirical 
correction of the curves from 3° trim to 0° trim may be seen 
from reference 2 to require a 10-percent reduction in the 
acceleration. An analysis of data obtained at the Langley 
impact. basin for impact at 0° trim indicates that correction 
of the time scale from 3° trim to 0° trim requires a 10-percent 
reduction in the time values. 

Hydrodynamic force.—The nodal-point acceleration ns 


obtained by the procedure discussed herein, represents the . 


hydrodynamic force, in multiples of the weight, applied to 
the flying boat. Since the experimental data do not provide 
measurement of the hydrodynamic force as such, direct com- 
parison of the theoretical force-curve results with experi- 
mental results is not permitted. Instead, & comparigon of 
the theoretical response of the structure with. the experi- 
mental response is necessary, and, if the agreement is ade- 


quate, it may be concluded that both the hydrodynamic and _ 


the structural actions are adequately represented. 

Before a study of the response of the structure is made, the 
theoretical effect of the response on the hydrodynamic force 
should be observed. This observation is made by comparing 
the force curve obtained for a mass ratio of 0.25 with the 
force curve for & rigid structure. Both curves are included 
in figure 11. The curve for the case of a rigid structure 


m ; 
(22-0) was obtained from reference 1 for 3° trim and cor- 
L 


rected to the conditions of the present example as previously 
indicated. Comparison of the curves for mass ratios of 0 
and 0.25 shows that in the present example the theoretical 
effect of the strvctural elasticity on the maximum hydro- 
dynamic force is to reduce it 15 percent. For the hypo- 
thetical condition of a concentrated wing mass located at a 
point in each semispan of a massless wing structure, the 
conditions of the present example would give theoretical 
reduction in the maximum hydrodynamic force due to 
structural elasticity of 44 percent. This result is indicated 
by comparison of the maximum of the curve in figure 11 for a 
mass ratio of 1.36 with the maximum of the curve for mass 
ratio of 0. 

The curves in figure 11 show only reduction of the hydro- 
dynamic force because of structural elasticity. . It should 
not be concluded, however, that. the effect is always in this 


direction. Figures 2 to 6 show that in some cases the hydro-' 


dynamic force is increased; the maximum increase which was 
calculated. was of the order of 12 percent. 

Hull acceleration. By combining the results for the two- 
mass case with approximation of the pitching action on the 
basis of & rigid structure, the following equation may be 
obtained for the hull acceleration at. different longitudinal 
stations: 

Ty =ni f +h (7) 
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= 40,000 pounds; fa 3.6 cycles per second. 


where 
Tr hull acceleration 
8 
f-ltg 
'8 distance from station to center of gravity 
k radius of gyration (12 ft in present example) 
L distance from resultant hydrodynamic force to 


center of gravity. 

The forebody length of the flying boat is 31.75 feet. Be- 
cause of bow effects, a length of 25 feet is assumed to_have a 
rectangular loading for this zero-trim impact and the result- 
ant force is located 12.5 feet forward of the step which leads 
to a value of / equal to 8 feet. . Use of the foregoing procedure 
to calculate time histories of the acceleration for two stations 
in the hull at which accelerometers were located gives the 
curves in figure 12. The maximum accelerations recorded at 
these stations are also shown; agreement with the computed 
maximum acceleration is good. The full experimental time 
history is not included because the film speed was not great 
enough to permit accurate determination of the shape of the 
time history. This factor, together with some uncertainty . 
in defining the exact instant of contact, prevents exact check 
of the time to reach maximum acceleration; thus, the experi- 
mental points in figure 12 are located at the theoretical 
time of maximum acceleration. 

Although the agreement of calculation with experiment in 
figure 12 is good, this agreement can be interpreted as con- 
firmation of the elastic action of the structure only to the | 
extent that disagreement of hull accelerations computed on 
the basis of & rigid structure can be shown.. Values of 
maximum accelerations computed for & rigid structure are 
included in figure 12. The disagreement with experiment 
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Figurer 12.— Variation of bow and step acceleration with time. 


is greater than for computations in which the elasticity of the 
structure is considered. The difference is small, however, 
because the effect of elasticity in increasing the hull accelera- 
tion for a given force is largely offset by the effect of elasticity 
in reducing the hydrodynamic force. Since the probable 
accuracy of available data must be considered, experimental 
proof of the theoretical effect of the elasticity on the hull 
acceleration cannot be claimed. Support for a conelusion in 
this respect, however, can be obtained by comparing the 
theoretical response of the wing with the experimental 
response. If the wing responds as assumed, the basie equa- 
tions require that the hydrodynamic force and hull accelera- 
tions be as calculated. 

Elastic axis.—Use of the nodal and oscillatory accelera- 
tions of the representative two-mass system to predict 
accelerations along the elastic axis of the wing requires 
consideration of the fact that the elastic axis of the wing is 
not at the center of gravity of the flying boat. An approxi- 
mate correction may be obtained by multiplying the results 
for the two-mass case by the factor f which is used in equation 
(7). This correction is not entirely consistent with that 
given by equation (7), but each approximation is considered 
more accurate for its particular case. An improvement to 


the. present correction which would change the results 2 to 3 - 
percent might be made, but the complication is not considered 


to be warranted. 

In the present example the value of f is 0.86. Application 
of this factor, of equation (B7) of appendix B, and of per- 
tinent structural data given in reference 3 gives the curves 
in figure 13 in representation of the acceleration time histories 
of the elastic axis of the wing for the hull, inboard-engine, 
~ nodal, outboard-engine, and tip stations. 

Wing torsion.— The torsion of the wing during impact may 
have substantial effect on the acceleration of engine and 
nacelle masses forward of the wimg. Use of the procedure 
which gave the acceleration time histories in figure 13 to 
calculate the acceleration at the engine gives a maximum 


~ deflection of a 


acceleration of 3g at the inboard engines and a maximum 
acceleration of 5.6g at the outboard engines. The relative 
magnitude of these values is in strong disagreement with 
the fact that during this impact partial failure of the inboard- 
engine mounts occurred, but the outboard-engine mounts 
were not damaged. 


Since the structural data used in the preceding computa-_ 


tion are for the actual fundamental mode, their use in impact 


. calculations involves assumption that the coupling between 


the torsion and bending corresponds to the coupling which 
exists in natural vibration. Actually, the torsional deflec-_ 
tions in impact are determined not merely by the bending 
deflections but also by the large nodal acceleration, whieh 
does not exist in natural vibration. 
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In the present example the results which have been caleu- ` 


lated for impact of a two-mass elastic system will be used to. 


predict the response of the engines. The procedure for 
doing this is to select the proper solution and then to adjust 


ihe acceleration and time scales of the two-mass solution so, 
that the maximum acceleration of the nodal point corre- 


sponds to the maximum acceleration of the elastic axis at the 
engine station. The acceleration of the upper mass of the 
two-mass system then represents the response of the engines; 
however, because of the eccentricity of the impact, an in- 
crement must be added. 


Data obtained from Guggenheim Aeronautical Labora-  _. — 
tory, California Institute of Technology include the torsional — = 


station inboard of the outboard engine for a 
given moment applied at the wing tip. For a flying boat 
the absence of cut-out for the landing gear tends to give a 
stiffer wing; therefore, in the present example, deflection 


measured on the wing of the landplane slightly inboard of | 


the outboard engine is considered applicable to the outboard- ~ 


engine station. Relative deflections between the inboard 
and outboard engines are estimated as follows: 

(1) Torsional deflection at inner engine equal to 1 unit 
due to each engine, or 2 units total 


(2) Torsional deflection between inboard and outboard 


engines equal to 2 units because of greater distance of flexure, 
boosted to 3 units because of increased flexibility of structure. 
On the basis of data from California Institute of Tech- 


nology the average of the static moments of the inboard and = 


outboard engines is taken as 22,000 foot-pounds. After 
determination of the torsional deflection at the outboard 
engine for this average moment applied at the wing tip, 


multiplication by the ratio 1/2 gives a value for the static 


deflection of the inboard engine and multiplication by the 
ratio 5/4 gives & value for the static deflection of the out- 
board engine. 

An approximation to the calculation of the response of the. 
engines to the total acceleration of the elastic axis is to 


neglect dynamic interaction between the engines and treat ~~ 


each engine as a single-mass oscillator having a natural fre- 


dure gives values of 7.6 cycles per second for the inboard 
engine and 5.1 cycles per second for the outboard engine. 
Use of these frequencies to determine t, and division of L . by 


quency determined by its static deflection. Such a proce- 
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FIGURE 18.— Variation of accelerations at stations along the wing elastio axis with time, 


t, (E equal in this case to the time to reach maximum acceler- 
ation of the curves in figure 13 for the station in question) 
leads to values of t,/t; for the inboard and outboard engines. 
The next step is to select a mass ratio which for these values 
of t,/t; has & shape of the nodal acceleration time history, 
which approximates the shape of the acceleration time history 
for the elastic axis at the station in question. In the present 
example the mass ratio of 0.25 is used. Seale factors for both 
the load and time scales are determined so that the maximum 
nodal acceleration for the two-mass solution will agree with 
the maximum acceleration and time to reach maximum 
acceleration of the elastic axis at the station in question. 
After these factors are applied to both the nodal and oscilla- 
tory curves for the selected mass ratio and time-period ratio, 
use of the results and equations (B7), (B5), and (B2) in 
appendix B to calculate the acceleration of the sprung mass 
of the two-mass system gives accelerations of tlie engines. 
Approximation and superposition of the pitching action on 


the basis of & rigid structure leads to the solid-line curves 
given in figure 14 for the accelerations of the engines. 

Also included in figure 14 is the design ultimate acceleration 
for the engine mounts. Comparison of the caleulated engine 
accelerations with this value shows agreement of the calcu- 
lation with the fact that partial failure of the inboard-engine 


. mounts occurred but tbe iiiad mounts were not 


damaged. 

During the impact an accelerometer was located at the 
outboard-engine station intermediate between the clastic 
axis and the engine. A calculated time history of the ueceler- 
ation at the accelerometer location is given in figure 14; this 
time history is based on linear interpolation between the 
computed accelerations at the engine and the clastic axis at. 
this station in accordance with the proportionate distances 
involved. The figure also includes the recorded maximum 
acceleration at this point and shows good agreement of the 
computed acceleration therewith. 
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FIGURE 14,— Variation of accelerations of engines with time. 


Acceleration time histories for the engines, computed on 
the basis of a rigid structure, are included in figure 14. The 
maximum accelerations computed on this basis do not agree 
with the structural failures which occurred. F eT 
the fact that the curves computed on the basis of a rigid 
structure reach a maximum at the same instant of time is in 
strong disagreement with experiment, which in this respect 
is in approximate agreement with the computations for an 
elastic structure. 

Difference between the ratio ¢,/&; for the inboard-engine 
and outboard-engine stations is primarily responsible for 
difference in the calculated response of the engines. Most 
of the difference in this ratio for the two stations is not due 
to difference in the sprung-engine frequency but is due to 
the greatly different time to reach maximum acceleration f; 
of the elastic axis. Agreement with experiment of the 
response calculated at these stations provides indirect con- 
firmation of the acceleration. time histories predicted for the 
elastic axis at these locations by the normal-mode method. 

For the impact, experimental data are not available for 
checking the tip acceleration, which is predicted on the 
basis of the normal-mode method, but the initial downward 
acceleration and the 12g maximum acceleration shown in 





Yt L3. 
a 


figure 13 for this station are in general agreement with, ] 


results recorded in severe impacts of other airplanes. Agree 


ment of the computed hull acceleration with experiment has 


already been shown; in an indirect manner all the curves in 
figure 13 exhibit satisfactory agreement with available ex- 
perimental data. Although the response of the engines is 


different from the response assumed in celeulating these ` 
curves, it appears that in practical use the two-mass solu- . 
tions given herein can be interpreted on the basis of the __ 
normal-mode method to obtain both the response of the 

` elastic axis and the hydrodynamic force. 


that in practical problems modification of these two results 


is not required when accelerations, of the engine different 


from the accelerations predicted by the normal-mode method 


Further, it appears ` 


are determined; the merit of this statement should be | 


independent of whether such modification is made by the 
method used herein or by another method. 

In the foregoing comparison agreement of calculation with 
experiment is obtained without consideration of the response 
of modes higher than the fundamental mode. If, when 
more data are available, it is shown that the response of the 
higher modes can be determined by treatment parallel to 


that given the fundamental mode, the two-mass solution 
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given herein can be used to predict their response by selecting 
a solution, for a mass ratio and natural frequency representa- 
tive of the higher mode, which has a nodal acceleration curve 
of T the same shape as the hydrodynamic-force 
curve 


CONCLUSIONS 


Theoretical solution of hydrodynamie impact of a hull 
mass connected by & spring to an upper mass, the results of 
calculations for wide ranges of mass ratio and natural fre- 
quency, and the use of these results in & comparison of 
theoretical data with test data for a Bight-test landing 
impact indicated that: —— : 

1. In. flying-boat impact the effect of the structural 
response on the hydrodynamic force might be substantial, 
the shape of the force time history might be considerably 
ehanged, and the maximum hy drodynamic force might be 
gei reduced or increased. 

2, The greatest reduction in hydrodynamic force occurred 
for the condition of large mass ratio and low value of spring 
constant. 


3. The normal-mode method was a practical means for 
determining the equivalent two-mass system which repre- 
sented the major elastic action of the airframe, for predicting 
the effect of this action on the hydrodynamic force, and for 
approximating accelerations along the elastic axis of the 
wing. 

4. The acceleration of engines contained in nacelles for- 
ward of the wing could not bé computed on the basis of 
coupling between torsion and bending as m the fundamental 
mode, but simple treatment of their response to the com- 
bined translation and oscillation of the clastic axis gave 
agreements with experimental accelerations and gave an 
explanation of partial failure of the inboard-engine mounts 
during impact. 


LANGLEY MEMORIAL AERONAUTICAL LABORATORY, 
NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS, 
LaNGLEY FrELD, Va., March 17, 1947. 
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APPENDIX A 


MATHEMATICAL EQUATIONS AND METHOD OF SOLUTION 


EQUATIONS OF MOTION 


The following equations of motion, which were derived 
from equation (30) in reference 2, are for fixed-trim impacts 
of a rigid prismatic float connected by a massless spring to 
2 rigid upper mass: 

Acceleration of lower mass normal to water surface in feet 
per second per second 





4a, d 
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Acceleration of sprung mass normal to water surface in 
feet per second per second 


(A1) 


" L. T — 

ys=7 Ayt me) Y H3 Ay (yo Kicosry] (A2) 
Acceleration of nodal point normal to water surface in feet 
per second per second 


- MryrtMsy 
pA abs (A3) 
The spring constant in pounds per foot of deflection can be 
expressed by the equation i 


< {rm mMms) fa " , 
Bv ML EMs e 
where 


sin?r 
Vn COS T 


T tn 7 
mdr (3571) [eds tan aL aada S ) 


jn, velocity of lower mass normal to water surface, ft/sec 
ye draft normal to water surface, ft 

8 angle of dead rise, radians 

r angle of trim, deg 

yz, initial velocity normal to water surface, ft/sec 
£z, initial velocity parallel to water surface, ft/sec 
mg sprung mass, slugs 

mr, lower mass, slugs 

f. natural bending frequency, cycles/sec 

p mass density of water, slugs/cu ft 


Kest T — 


® ee ®® © 6 ® ® © © © © 0000000 


COMPUTING DIRECTIONS 


A sample data sheet is given as table I. In this table de. 
numbers in circles refer to rows; the circled numbers under m 
the row headings refer to computed values to be used for 
the computations. In the first column the time is equal to 
zero, in the second column the time is AZ, and in successive | 
columns the time is 2 Af, 3 At, and so forth. Each row is 
computed in sequence for any given column before any row 


is computed for the next column, except for the first column. 


TABLE I.—DATA SHEET—GENERAL TERMS 
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In the first column the assumed value of +: is the given value 
of Yro 
@), €3, and @3 which in this case are also yz, Lower-case p 
as a subscript on a row number refers to the computed value 
in the indicated row of the preceding column. 

Each row heading indicates the general operation to be 
performed. When the data sheet is set up, numerical values 
should be substituted in the row headings for the expressions 
that are constant for a specific case. The constants for the 
sample computation, given in table II, are as follows: 


R EE ED eae NONSE oe ieee MC SHE us. 107, 664 
Pre ee i et Ra oie Ec UM E eo 5 E A 
S. EEEN LE cetera eee naa oa ee uidi sie oan a -183, 919 
B, SEE - - 22-5 
T, dep. coscienza dora dial Uus eee aes tee m = Tk a8 
UNE G1Ó —— — P——— — 30. 6673 
JL —— aacr anc a BÀ 1589 
[P MEME ce KNEE 715. 217 
We evasit. mE ET edente uam cuenssmuuui 525. 776 
fs, cycles/sec__.-.-_.~- nx ee OE ene a 
p, slugs/cu ft___..-..--.~ PP Sede dix Bee coe eel ee, S [098 
Al, S00 iz ced RE sine cuneate E b EE DOS 





For most solutions a time increment At of 0.005 second is 
satisfactory. (The value of At may be varied by considering 


All other values are zero except the values for rows . 


the time for a given mass to reach a maximum acceleration.) 

' The number of significant figures to be used should be 
chosen on the basis of the computing equipment available 
and the accuracy desired. Comparison of a solution com- 
puted with four significant figures with results which had 
been obtained with six significant figures gave a difference of 
about 1 percent in the maximum acceleration. 

As the computations proceed, the lower-mass acceleration 
(row @), the sprung-mass acceleration (row @), and the 
nodal-point acceleration (row @@), all in g units, should be 
plotted against time.(row (D) in seconds! Each time & 
column is completed the new points should be added to the 
plot. This plot is the only brief method of checking on the 
accuracy of the. computation until sufficient solutions are 
obtained to permit cross-plotting. If the points do not lie 
on a smooth curve, an error has been made in the computa- 
tions. If an error is made, every value computed thereafter 
contains the error and therefore great care must be taken. 

The number of columns required for the computations 
depends on the number of columns required for rows @ and 
62 to pass their respective maximums. At least three 
columns should be computed beyond the column in which 
row @ reaches its maximum value. 


TABLE IL.—SAMPLE COMPUTATION 
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Time At, sec (0.005)... i es 0. 005 . 0. 0. 015 0. i . 08 i 
Be aa O 20. 6673 20. 6673 20. 6420 20. 5287 20. 2688 19. 8200 19. 1964 18. 3859 
(4), 0.005) +@,_----------- . O | .Q, 103337 |. Q. 206579 | 0.309380.) 0. 411128 | 0.511025 | 0.608174 | 0. 701700 
STO aaa ——1... -. | 0.010679 | 0.042675 | 0.095716 | 0.169028 | 0. 2611147 | .0. 360876 | 0. 192383 
IMS PRU IPEA O | 0.001104 |. 0. 008816 | 0.020613 | 0.069491 | 0.133453 | 0.224949 | 0. 345505 
TT T . O | 0.000114 | 0.001821 | 0.008102 | 0.028570 | 0.008108 | 0.136808 | 0. 242441 
(9 1:6 297306) 9 X 401. 251... 0 2601.07 |  10612.5 23586. 5 40782. 0 60767. 4 81002. 9 | 101248. 0 
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MOMMSEN -O | 0.003105 | 0.031778 | 0.140730 | 0.422190 | 0. 997595 2. 01448 3. 63288 
Bren Tore 0| 0.000027 | 0.000160 | 0.000506} 0.001168 | 0.002213 | 0.003821 | 0.005977 
( 2. 00281. S ..0 0. . . O | 0.000002 | 0.000006 | 0.000015 | 0.000030 | 0.000054 
ig m "EU cM O | 0.014427 | 0.043282 | 0.108204 | 0.216409 | 0. 380525 
HEN .  410| @ 007222 |—0.023903 [—0. 319925 |—0. 138058 | —3. 63246 | —7. 88831 | —14 7856 
PA ERE TO 0] 2662.16.) 106089.) 235383 40581. 7 60220. 6 80415. 5 90028. 0 
133.919 X ®©) +525.776__._._. -0 525. 024 526. 957,| 529.742 585. 082 543. 648 555. 901 572. 046 
VORNE ..0 | —5. 06187 | —20. 1324 | —44. 4335 | —75. 8420 | — 110.771 | —144. 658 | —173. 112 
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849... aar hao 20. 6873 20. 6420 20. 5287 20. 2688 | '19. 8290 19. 1964 18.3859 | 17.4355 
—0. 001398(@— 48)]/32.2.. ~- ..0.| 0.000047 |—0. 000156 |—0. 002093 |—0. 008696 |—0. 023740 |—0. 051552 |—0. 096297 
322. 0 |—0.157201 |—0.625230 | —1. 37992 | —2 35534 | —8. 44009 | —4. 49218 | —5. 37615 
(625.776 ET B) |. . 0 |~0. 066575 |—0.264083 |—0. 585844 | —1.00291 | —1. 47116 | —1.93300.| —2 33323 
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APPENDIX B 


RELATIONS BETWEEN TWO-MASS SYSTEM AND REPRESENTED STRUCTURAL MODE 


The sum of the masses m; and ms shown in figure 1 must 
equal the gross mass of the represented airplane in order to 
obtain the proper nodal acceleration. For the hypothetical 
limit condition in which the wing mass is concentrated at a 
single point in. each semispan, m; is the actual hull mass and 
mg is the actual wing mass. In order to take into account 
the more complex nature of the structural action for a par- 
ticular mode, the determination of the ratio of these masses 
is necessary so that the vibrational energy of the simplified 
and represented systems are equal for the same vibrational 
amplitudes of mj and the actual nu or point of force 
applieation. 

On the basis of the theory of saias, for which equa- 
tions are included in reference 3, the vibration energy E of 
the two-mass system is given by the equation 

el mier Msgs (B1) 
where 


vy, Vibrational amplitude of mz relative to nodal point of 
system 

es vibrational amplitude of mg relative to nodal point of 
system 

e natural frequency 


Since for the two-mass system the node is at the center of 
gravity, 
(B2) 


mMrYr—=—MNsPs 


If equation (B1) and equation (B2) are eombined, and since 
the total mass m is equal to the sum of the masses mg and 
m;, the following equation can be obtained for the vibrational 
energy E of the two-mass system: 

E 1 mM r, 2 


— c 
wo” 2 ma Y 


(B3) 





The vibrational energy of the represented mode is a function 
of the spanwise mass distribution and mode shape. On the 
basis that the semispan of the airplane is divided into 7 sec- 
tions or stations. the vibrational energy of the mode can be 
written as follows: 


E 


1 
wv 2 (miei tme +t meo t ... +m) 


where 


m, mass of jth spanwise section; value for semispan doubled 
to represent the entire span 

€; deflection of mass at 7th spanwise section relative to 
nodal point (A, ta) 

- hy deflection of elastic axis at jth spanwise deflection rela- 
tiye to nodal point 

zx chordwise distance from elastic axis to effective mass 
center 

æ; torsional deflection at 7th spanwise station 


(B4) 


Equality of the vibrational amplitude of the lower mass of 
the simplified system to the vibrational amplitude of the 
hull or fuselage of the flying boat or airplane relative to the 
nodal point of the represented mode is expressed by the 


equation 
EL = Pk 


where 


g& deflection of hull or fuselage of fying boat or. airplane 
relative to nodal point of represented mode 


The requirement of equal energy of the simplified and 
represented systems for the condition expressed by equation 


(B5) gives combination of equations (B3), (B4), and (B5) 


to obtain the following equation for the mass ratio of the 
two-mass system: ' 


ms _ me 
mr " mig F Maga? mio! + ee 


E ne fe 6 
Tm, n) 

Computation of the natural frequency of wing modes has 
received a great deal of attention in connection with study - 
of wing flutter and need not be treated herein. Incidental 
to calculation of the natural frequency, & mode shape is 


attained which, together with knowledge of the mass distri- 


bution, permits use of equation (B6). In cases in which the 
wing has been constructed, the mode shape and natural 
frequency may be determined experimentally. Equation 
(A4) in appendix A of the present report permits computa- 
tion of the spring constant which for & given mass ratio of 
the simplified system gives the required natural frequency. 
After the accelerations of the two masses of the simplified 
system have been computed, equation (À3) of appendix A 
fixes the magnitude of the nodal acceleration. The diference 
between the nodal acceleration and the hull, fuselage, or 
float acceleration can be taken 85 a measure of the oscillatory 
acceleration. On the basis that the structure deflects in the 
mode used in determining the equivalent two-mass system, 
the acceleration at any point is given by the equation 





n,—ncbneb = -.- - (BT) 
eh 


The foregoing equations, with consideration of rotatory 
inertia and energy, may be applied to the case of a tip float 
attached to a flexible wing if the stiffness and mass distribu- 
tion are known and if a manner of structural deflection is 


assumed. 
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